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DYNAMIC ANALYSIS OF EVOLUTIVE CONSERVATIVE 
SYSTEMS. DISCUSSION OF "EIGENMODE 
CROSSINGS" 


H . , J. -P . Morand 


ANNOTATION 


Following an analysis of the close link that » 

exists between the symmetries of a dynamic system and jj 

the multiplicity of its natural vibration fre- jj 

quencies, we show , on variational grounds, that l 

for a system with a fixed symmetry, the natural 
frequencies associated with natural modes of a 
g: ven type of symmetry do not cross each other 

during the evolution of this system. , 

The theoretical method has been put to use 
numerically and applied to the analysis of the 
evolution of the axisymmetric hydroelastic modes 
of the "Ariane" launcher during first-stage 
combustion. 
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DYNAMIC ANALYSIS OP EVOLUTIVE CONSERVATIVE 
SYSTEMS*. DISCUSSION OF "EIGENMODE 
CROSSINGS" 


H., J.-P. Morand** 


INTRODUCTION 


To set the stage, let us consider the evolution of the /2* ** 
longitudinal modes of vibration of a liquid propellant 
launcher during combustion of one stage. In practice, it 
is convenient to characterize the conf iquration of such an 
evolutive dynamic system by means of the draining rate of 
the propellants in the functioning stage . We might think 
that the change, a priori considerable, of the natural 
frequencies resulting from changes in the mass of the vehicle 
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is necessarily accompanied by frequencies crossing each other 
(Figure 1). However, when the experimenter's track the natural 
frequencies of an evolutive system, they observe, in fact, 
events of the type depicted in Figure 2 and that we have 
agreed to call "modal interaction". 


frequency 




Fig. 1 - Frequency Crossing. Fig. 2 


"Modal interaction". 


One of the goals of this study is to show -as if there was 
a need!- that experimental reality matches theoretical reality 
and that, with the exception of relatively academic cases 
that we will discuss, the eigenmode crossings must be 
considered as "accidental" with respect to modal interaction 
events . 


If we were only discussing the physical reality of eigenmode 
crossings, our intent could appear to be somewhat academic 
and provide only a principle for plotting curves describing 
the evolution of frequencies. However, it happens that during 
modal interaction, the ability to excite a mode of vibration 
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with a given type of excitation can undergo "catastrophic" 

variations: And this excitability measured by the effective 

mass u (defined below) proves to be playing an essential 
© 

function in predicting pogo instabilities of launchers using 
liquid propellant for which we now recall the hydraulic and 
mechanical origins (ref. 1). 

In flight, the thrust provided by the engines can fluctuate 
slightly on either side jf an average value. These 
fluctuations induce vibrations in the launcher structure, 
thus subjecting it to accelerations around natural resonant 
frequencies. These accelerations of the structure are in 
turn coupled to the system supplying fuels which oscillate 
at the resonant frequencies of the tanks and the supply lines 
to the engines. These last oscillations generate pressure 
peaks in the propellants. These pressure peaks act on the 
engines by causing the propellant flows to vary, thus causing 
thrust fluctuations at the frequency of the pressure 
fluctuations. If the natural frequencies of the sLructure 
and of the propulsion systems coincide, the pressure 
fluctuation resulting from the thrust fluctuation increases. 
It can even diverge if the thrust fluctuation generates a 
pressure fluctuation greater than the one that created it: it 
is the pogo instability. 

This basic analysis clearly shows why the investigation of 
the pogo stability loop requires, in particular, that the 
frequencies and the natural forms of hydroelastic vibrations 
of the launcher be known. The parameter measuring the 


sensitivity of a vibration mode of the tank system to an 

excitation at the anchor point M for the engine structures 

is the "effective mass" u defined by: 

g u ~ 

* » 2 m ' ( i ) 

with Ug designating the generalized mass of the vibration 

2 

mode for an arbitrary normalization and u (M) being the square 
of the modal amplitude of the displacement of the anchor point 

M. 

Now, we are going to analyze the interactions of modes, as 
well as resulting effective mass variations, in the general 
case of an evolutive conservative dynamic system. By 
conservative dynamic system, we essentially mean a system 
with harmonic variations described by a variational principle 
involving a potential energy and a kinetic energy that can 
be expressed as quadratic forms of acceptable displacements. 

N. B. The qualifier "conservative" implies that we neglect 
any dissipating effect in the vibration study and that, in 
the case of a system varying as a function of time, the 
mechanical characteristics of the system (stiffness, mass) 
can be considered as constant over several periods of the 
vibrations under study. In concrete terms, in the case of 
the Ariane launcher, the combustion time for the first stage 
is 140 seconds whereas the period of the vibrations studied 
is less than one tenth of a second; it is therefore proper 
in such a situation to compute the vibrations at a given time 
by freezing the draining rate t which describes the 


corresponding instantaneous configuration. 


During the first part of this study we will try to define, 
as generally as possible, the forms and natural frequencies 
of such a dynamic system. 

We will consider the close link that exists between the 
multiplicity of the natural frequencies and the symmetries 
of a dynamic system. This indispensable tendency to decrease 
will enable us, on the one hand, to show theoretical cases 
of natural frequency crossings and, on the other, to state 
in a constructive manner the problem of the evolution of the 
vibration modes of a system. 

We will then show how a so-called "two-mode" approach to the 
solution of the variational problem governing the natural 
forms and frequencies makes it possible to prevent the 
expected frequency crossings in the one-mode approach. 

We will also give a geometric interpretation of the 
interaction of the two modes and of the remarkable nature 
of the effective mass that can result from it. 

During the second part of this report, we will illustrate 
the general analysis outlined above by applying it to the 
discussion of the evolution of the hydroelastic modes of the 
Ariane launcher during first-stage combustion. 



We will present some numerical results involving an analytical 
simplification of the evolution problem of the vibration modes 
of a liquid launcher valid for the investigation of the 
configurations close to a configuration for which we know 
the modes of vibration (in the present case these natural 
modes have been computed using the method of finite elements). 


1. - VARIATIONAL MODELING OF AN EVOLUTIVE DYNAMIC SYSTEM 

1.1. Varational formulation of a problem of harmonic 
vibrations . 

In order to introduce the general principles that enable us 
to describe a general conservative dynamic system, let us 
consider, as an example, the classical problem of the 
determination of the elastic vibration modes of a body (ref. 
2) occupying a domain n (volume) and embedded over a portion 
r of its boundary 0fl. 

By definition, a natural mode of vibration is a stationary 
solution of the problem of small motions of the body about 
its static equilibrium position. If U(M,.t) designates the 
displacement at time £ of the particle of matter located at 
point M at equilibrium, we become interested in the solutions 
of the dynamic problem of the form: 

*11 (M, 0 = u(M) cos (c of + </>) ( 2 ) 

with <p independent of M. Such a solution, when it exists, 
expresses a harmonic vibration, in phase, of every point of 
the elastic body. 


For a body with finite dimensions we can show that such 
solutions exist for certain discrete values of w: u>o»Wi,»Wa , 
... tu^..., that constitute a series of numbers . w n is called 
a natural pulse and the associated field u n ( M) is called a 
natural form. 

In the particular problem considered, if we designate with 
u q ( a=l ,2,3), the cartesian components of a certain natural 
form u. , the cartesian components of the Cauchy stress 

tensor considered as a function of the displacements by means 
of the behavior equation for the hyperelastic material, 

X the cartesian coordinates of a point in space, and n the 

Ot 

components of the unit line perpendicular to an and external 
to n , w the natural pulse, the couple (it,w) satisfies the 
equations : 


+ 00J 2 U =o at every point of 11(3) 

3 X, 

= o over an - r (4) 

u c* = 0 over r (a=l,2,3)) (5) 

(p = density per unit volume of the material) 

It is convenient to write the preceding problem in an 
equivalent form by means of the virtual works principle. 

For this, we introduce the class C of the fields of regular /4 
vectors defined over n and equal to zero over r. The dynamic 
problem is now stated as follows: 

Find w and ueC so as to satisfy the equation 


r ? 
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regardless of the 6u€C 


f o afi (u)e aj} (8u) - w 2 / ' p u a (8u) a - 0 

bu a 8Un \ 

(with e afi (u) = — + - — ) 
p dx 0 dx a ' 


(6) 


We say that C is the class of the "admissible forms" of the 
problem and that the condition at the limits (5) satisfied 
by the u elements of this class is essential. 


As for the conditions at the limits (4), we observe that it 
no longer appears explicitely in equation (6): indeed, we 

can show that it results from satisfying the variational 
principle (6). This condition at the limits (4) which is 
included in the variational principle (6) is qualified as 
natural . 

2 

Finally, let us consider the ratio [ id ] 


[co 2 l = 


1 

2 



( 7 ) 


considered as a function of „ <= S; \ o a(3 e 0a and ^ f pu 2 

constitute quadratic forms of u£C which have a direct 

b 

physical meaning. Then, let us introduce into expression 
2 

(7) for a) a field of the type U- + €6tt, being a solution 
of (6) for a "natural value \ n = w 2 w " and 6 u any class C 
field. We can then check, by noting that 

f a ct(s(“) a (8u) anc f P« 5u are syntnetrical bilinear forms, that the term 
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in e of the development of Ew 2 ] is zero. Stated otherwise, 
we note that; the function (7) defined over & is stationary 
at points u n of C and that the value taken by this function 
at such points is equal to The reciprocal is direct; to 
give a general problem statement that we will use as a 
starting point for the analysis of the conservative dynamic 
systems, we introduce two symmetrical bilinear forms 
K[u,v) &Nl(u,v) defined for a class of acceptable forms C 
verifying certain conditions at the essential limits. We 
will call the natural form u~ n of the dynamic system 
described by K and M any stationary point of the function 
[to 2 ] defined below: 


, 2i K(u , u) 

' M(u,u) 


6[co 2 )-0 


6 S 


(8 


(formal way of writing the condition for the standing state) 

The value of [w 2 ] at such points of C is the square of the 
associated natural pulse to . (For the general study of the 
dynamic systems, we will use the standard printing character 
u. to designate an acceptable form of C ) . 

1.2. Classification of the types of an evolution problem. 

We will now consider a conservative system depending 
continuously on a real parameter t and we propose to analyze 
how the natural forms and the natural pulses evolve as 
a function of t. 


We are led to distinguish between two main categories of 
mathematical problems. 

Category A: "evolution of operators". 

By definition, it includes cases where the 0 class of the 
admissible forms is independent of t and where operators K 
and M are dependent on t . 

Let us note that stating that C is constant means in practice 
that the domains (volume, surface, , . , , )over which the unknown 
fields of the problem generally noted as "a" are defined are 
constant, and also that the conditions at the essential limits 
satisfied by a are constant. 

As a typical example of a problem of this category, we can 
mention the problem of the determination of the vibration 
modes of a homogeneous body occupying a volume fi and embedded 
over a fixed portion r of its boundary as a function of the 
density of this body, the elastic characteristics (behavior 
equation) being assumed as constant. Moreover, in this case, 
only the operator M(u,v)~ f n P uv depends (linearly) on 

t = p. 

Category B: "evolution of the class C". 

We group in this category the problems of the perturbation 
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of conditions at the essential limits (the domains remaining 


constant ) and the problems of the perturbation of the 

boundaries of domains for which the domains depend 

continuously on a parameter t. 

In the rest of this report, we will only consider problems 
of the first category. Yet, the example of liquid-propellant 
launchers that we propose to treat seems to typically belong 
to the second category! Nevertheless, we will show that 
through a judicious choice of the unknown fields, we can go 
back to a problem of the evolution of operators. 

The first question that we wish to answer is the following: 
when x varies, can two distinct natural frequencies for 
t= 0 that remain continuous as a function of x merge for a 
particular value of x ? In other words, can there be a 
frequency crossing in the evolution curves for the natural 
frequencies of the dynamic system as a function of x ? 

We observe that such a crossing (Figure 1) must be accompanied 
by a degeneracy (or multiplicity) of a natural frequency. 

t 

Therefore, we are naturally led to analyze under which 
conditions such a degeneracy is predictable. 

1.3. Degeneracy and symmetry 

We can say that as a general rule the degeneracy of the 
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natural frequencies of a system is closely linked to the 
symmetries of this system. To define this link, we will now 
recall a few results from the "group theory" (ref. 3). 

Let G be the group of the symmetries that keeps a dynamic 
system constant. In concrete terms, in tne case of a bounded 
system, the possible symmetries are by necessity one of the 
types: rotations about an axis, planes of symmetry or center 
of symmetry, with the exception of translations for obvious 
reasons . 

Let us then designate as E^ the entire set of natural forms 
corresponding to a natural pulse w. Let us choose within 
E a base cfo, <\> 3 , . , . , ^ of natural forms so that any natural 
form -tj)€;c w can be represented in a unique manner through a 
linear combination of these n forms (n is by definition the 
multiplicity of w ) . 

Under these conditions, to state that the system studied is 
fixed with a a € G symmetry means that if ip is a natural form, 
the o <|> form obtained "by forcing the symmetry operation a on 
the natural form " is also a natural form for the same 
frequency and, consequently, belongs to E^ . For any^eE^. 
We can therefore write a4» in the form of a linear combination 

Of <j> 2 > tylf • • • 

We thus verify that for any element of symmetry a^G there 
is a certain linear substitution denoted as To within E 


Such an identity a : T» To which associates any oeleme«t of 
the G group, a linear substitution in the vector space E^is 
called a linear representation of the G group within E^. 

We can then attempt to decompose E^ into a sum of separate 
spaces E^^ f E^^ + . . . v-E^ ^ which are independently constant 
within G. 


We show that there is such a unique ultimate decomposition 
and we say that it provides the reduction of the 

representation T inside E, 

The identity a : T**Tj induced in each "irreducible sub-space" 

_ ( J ) 

is called an irreducible representation. 

_ ( J ) 

The dimension of the c' J ' space is called the dimension of 

co ^ 

the corresponding irreducible representation. 

The theory of the linear representation of the groups has 
for framework the vector spaces on the body of the C 

complexes, and to an irreducible representation of dimension 
4 by C there can correspond some physically irreducible 
representations of dimension 4 or 24 in the vector spaces 
on the body of the real numbers like the E^ space considered 
above . 

This remark is significant when this theory is to be applied 


to the classification of the natural modes of a conservative 


system and to the anticipation of the multiplicity of the 
natural frequencies. 

In the discussion that will follow we will use the following 
two results: 

The natural modes of a conservative dynamic system depicting 
a group with G symmetry are classified by type depending on 
the irreducible representations /yin 

The ft degeneracy of a natural space E of type T • is simply 

(i) j 

connected to the dimension 4; of the irrreducible 

J 

representation considered (ft = 4- or ft =24.- depending on the 

J j 

cases ) . 

In order to illustrate the purely geometrical meaning of the 
normal degeneracies of the natural frequencies, let us mention 
two theorems valid for the finite groups. 

- Burnside's Theorem: the sum of the squares of the dimensions 
of all the distinct irreducible representations of a finite 
group is equal to the number of elements of that group. 

The number of irreducible representations of a finite group 
is equal to the number of classes of that group (let us recall 
that we call a class from a G group a sub-set of G composed /]_ 
of elements that are "conjugate" elements of a given element 
a, that is to say generated by x. ^ ox. x. covering G; the 
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set of classes constitutes a partition of the group) » 

The preceding discussion enables us to consider some cases 
of evolution of the dynamic systems necessarily leading to 
theoretical mode crossings: that is the case of systems such 
that for a certain value of x , a sudden change of symmetry 
takes place. 

For example, let us consider the modes of vibration of a 
rectangular plate with sides a and 6 as a function of - = x 
(a constant). We observe that for x = 1 (square plate) the 
symmetry group of the structure acquires a fourth-order axis 
of rotation perpendicular to the plate at its center. But 
it happens that the symmetry group of the square plate has 
real irreducible two-dimensional representations whereas the 
symmetry group of the rectangular plate allows only real 
irreducible one-dimensional representations. Consequently, 
we expect, on the basis of these purely geometrical 
considerations, that for x - 1 , certain vibration frequencies 
of the rectangular plate merge in pairs. 

- We eliminate in the rest of our report such a case of 
theoretical crossings by considering only evolutive dynamic 
systems that have a symmetry group independent of x. 

In accordance with a remark made above, we can classify the 
various natural modes of a dynamic system according to the 
various irreducible representations of the symmetry group. 


- We can demonstrate that the determination of the natural 
modes of distinct types T* can take place by means of distinct 
variational formulas of type (8); in addition, we can arrange 
it so that each natural frequency with a given natural 
multiplicity (which only depends on the irreducible 

representation T ■ ) appears only once among the solutions of 
this variational formula. 

We will not consider here the considerable practical 
consequences of these last statements with regard to the 
computation of the vibration (or buckling) modes of structures 
displaying symmetries (repetitive, for example); we will 
limit ourselves to draw the following conclusion that is 
useful for the purpose of this report: 

Given an evolutive dynamic system with a constant geometry 

6; then, the natural modes and frequencies of distinct types 

Ty being governed by distinct mathematical problems, nothing 

prevents the natural frequencies of type T< to "cross" natural 

j 

frequencies of type Tj. 

1.4. Definition of the evolution problem 

In order to set aside from the discussion the two cases of 
theoretical crossings analyzed above, when we again talk of 
the evolution of natural frequencies within this report, we 
will mandatorily imply that we are considering the evolution 



of the natural frequencies and forms of a given T* type for 
a dynamic system of a fixed symmetry group G . 

Let us consider a diagram showing the evolution of natural 
frequencies as a function of a parameter t as depicted in 
Figure 3 . 


We propose to analyze the evolution of the natural modes and 
frequencies in the neighborhood of a value of x that we choose 
as equal to zero, thanks to a judicious translation of the 
parameter x . 


Figure 3 . - Evolution 
of natural frequencies 


r- ; 



frequency 



evolution parameter 


In accordance with a remark in the preceding paragraph, we 
will assume, without the restriction of a generality, that 
formula (8) has been selected so that the natural frequencies 
are not degenerated (except perhaps at frequency crossings, 
if they exist I ? ) . 


In order to simplify notations, we will also assume that 


only the operator M is dependent of t and we will even limit 
ourselves to a linear dependency (these two simplifications 
are optional and the reader will be able to adapt the 
following analysis to the general case). We will also use 
the notation X . for the natural value that is the square of 
a natural pulse w . . 

'C 

For t =0 the natural values X ^ are governed, as well as the 

* 

. ^ I 1 

corresponding varational forms u, ,by the variational 

'Ci 

principle : 


0 = 5 [XJ = 5 


Kju , u) 
M(u,u) ; 


u e s 


(9) 


For t ^ 0 the natural values X* are governed, as well as the 

'C » 

new forms obtained u.* f by the variational principle 


0 = 5 [X 141 ] a 5 


K{ii , u) , 

M(u , ii) + t V{u , u) ’ 


wee 


( 10 ) 


i\) is the finite perturbation constituted by the increase 
of the quadratic form M(u,u).- 


We assume that for x=0, problem (9) is resolved. It is 
convenient for the considerations that will follow to assume 
that the ^ series is normalized so that M =1 ; 
indeed, if ^ is a non-normalized series, it is sufficient 
to define ^ / (M ^ M )^ ^ . in consideration 
of this convention, the orthogonality relations through the 
solutions of (9) are written 



(6 •' Knonecker's symbol) 
h-j 
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( 11 ) 


1.5. Variational resolution of the perturbed problem 
We wish to resolve problem (10). We then assume that the 
series makes it possible to represent with as much 
accuracy as we want an admissible form it & C . /9, 

Therefore we can look for the solution of problem (10) in 
the form of a linear combination 

E (12) 

i-i 

the a- values being determined by -t-he condition of 

A# 

stationarity of [X*]. It is proper to verify that the "best 
values of eg 11 for a given N satisfy the equation at the 
following natural values where X* is the corresponding natural 
value . 


S„ a, - X* [5„ + t Fo 4°> . uf)] tt, = o (13 ) 

(• i and j varying from 1 to N). 

The preceding method of resolution makes it possible, through 
the introduction of an increasing number N of natural forms 
of the uj^ series, to determine with an arbitrary accuracy 


the solution of the perturbed dynamic problem. 
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Also, let us note that [X*] being stationary at the u . * points 
of C , any error of order e in the natural form u* is not 
accompanied by an error in e z in the natural value X* , 

A* 

Therefore, it is reasonable to expect accurate estimates of 
the natural frequencies, even with relatively rough 
approximations of the natural forms (ref. 4). 

1.6. One-mode approximation 

The simplest approximation that we could make to solve (10) 
consists of stating that ul^is, as a first approximation, 
a natural form of the perturbed system (10). The correspond- 
ing natural value is then obtained by injecting the "test 
function" into the expression for the Rayleigh quotient 

C X*]. This estimate will be annotated x*^J taking (11) 
into account, we therefore have: 

>o) ^ 0) 

" 1 +rK(«®,«j°>) < 14 > 

By repeating this computation for each natural form u!^, 
we obtain for each value of T what we will call the estimate 
of the evaluation of the natural values in the one-mode 
approximation. In Figure 4, we have depicted a typical 
variation of these estimates; let us assume that the 
examination of the approximate evolution curves thus obtained 
leads to a frequency crossing for t = x 0 (Figure 4). 

In the neighborhood of such a value, it is legitimate to 
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question the validity of the one-mode approximation. We are 
then led to Improve our method for estimating the solutions 
of equation (10) as follows; 


1.7. Two-mode approximation 


r 

05 - / 


We wish to examine "under a magnifying glass" the area of 
the diagram (Figure 4) enclosed by the dotted line. The 
simplest refinement of the preceding model consists of looking 


for the solution of equation 
i and j being the modes involved 



Figure 4. -Evolution of the 
natural values in the "one- 
mode" approximation. 


(10) in the form a>U; + a- j , 
by this crossing possibility. / 10 



Figure 5. -Evolution of 
the' natural values in 
the" two-mode" approxi- 
mation . 


The best a ; and a- verify equation (13) which is explicitly 
x. j 

written : 


23 



> + T V,, T Vy 



0 


k (0) 


“/ J 


= X 


r v a 1 + r v i, 



hi 




(15 


(we have used the abbreviated notation V • • for 

-Cj 


V[u[°\u { j 


0 J 


l) 


We note immediately that the preceding couple of equations 
is formally identical to the equations governing the 
vibrations of a system composed of two coupled harmonic 
oscillators. The natural frequencies of the two oscillators 
in question are precisely the frequencies obtained in the 
one-mode approximation according to (14). 


The detailed discussion of the problem at the natural values 
being elementary, we will only indicate the resulting 
consequences within the scope of our analysis; they are easily 
understood thanks to the analogy that have just mentioned. 


As long as the natural frequencies of the oscillators are 
sufficiently removed, the effect of the xt /. • coupling on the 
natural modes and frequencies is negligible: in our problem, 
that means that outside the area bounded by the dotted line, 
the estimates of the frequencies in the one-mode approach 
are very close to the exact values. 


On the other hand, the xl / - * interaction plays a crucial role 
"when the frequencies of the two harmonic oscillators are 
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very close": In fact, the Interaction effect is to reveal 

a "forbidden frequency range" (or "gap" in English-language 
terminology) where the one-mode approach allowed for a 
prediction of a frequency crossing 1 

We should now check that the simultaneously consideration 

of an increasing number of modes to solve (10) does 

not modify the qualitative conclusions depicted in Figure 
5. There also remains the analysis of the interaction / II 

situations with more than one mode but this is outside the 

scope of the work that we have set for ourselves here. 

The entire preceding discussion bears on the fact that when 

t(/ [a ® ^ reaches zero, it cannot do so accidentally. 

~ j 

If we had not taken care to consider only modes of a T > type 

3 

of a system of symmetry G , it would have been otherwise. 
Indeed, V is by hypothesis an operator that is constant within 
G (since it respects the fixed symmetries of the system 
considered); but there is a theorem in theoretical physics under 
the name of "selection rule" which is stated as follows: 

Given two natural forms u^ aj and belonging to two 

natural sub-spaces that transform into two different 
irreducible representations T a and of the G group, and 
given a bilinear operator that is constant within G ; under 
these conditions 1/ [u^ a la L ^ ) is equal to zero. 


It is from this property that we have the uncoupling of the 
mathematical problems caused by the search for the modes 
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of distinct types that we have previoulsy mentioned, Let 

us return to the physical interpretation of Figure 5. Xf 

we "follow 11 the natural form u*(t), we note that outside the 

i 

area defined by the dotted line, this form differs very little 
from tty(O) as long as x < t 0 and that it is very close to 
u^(0) for t > t 0 . On the other hand, for x - t 0 it has the 

hybrid characteristics of u>(Q) and a*(0). 

The reversal of the physical characteristics of the vibration 
modes takes place naturally for the solution corresponding 
to the lower branch 

This reversal of physical characteristics of the vibration 
modes of a structure within an interval of variation of the 

parameter t which defines its configuration is a 

characteristic of a two-mode interaction. 

1.8. Geometrical Interpretation 

We can give the following geometrical interpretation of an 
interaction event which will be useful in discussing the 
evolution of the effective mass p . When using the approach 

consisting of having the dependency in t of the dynamic system 

rely entirely on M( u,u), it is convenient to consider that 
the u«(t } series constitutes a base of C that is orthogonal 
to the direction of the scalar product K( u.,v ) (on the condition 
of removing from the discussion the rigid body modes that 
have a zero "length" for the "standard" associated with this 



scalar product). 


We can then condense the description of the evolution of the 

two natural mode solutions of (15) by saying that during the 

interaction, the directions of u#( t )u? ( t ) perform a 

i 

simultaneous rotation of about 90 degrees while remaining 
perpendicular to each other and while remaining located 
within the plane defined by u^(0) and tty ( 0). 

1,9. Evolution of the effective mass during a modal interaction. 

We wish to show that during a two-mode interaction, the 

effective mass can become infinite for one of the modes. 

Let us return to the simple case of a structure excited at 

a point M along a given direction A, which was useful to us 

in introducing the concept of an effective mass (1). Let 

us then consider within C all of the admissible forms u. such 

that the displacement a(M) along the direction A of the 

excitation force is equal to zero. All of these represent 

a hyperplane of C . If we represent the plotted line V of /12 

this hyperplane in the plane generated by u. ( 0) and u»(0) 

x. J 

(Figure 6) we note that during the rotation as a unit of the 
couple of the directions of IL A t) u < ( x ) , the end of one of 
the two representative vectors crosses the line V for a given 
value of t . Consequently, for this value of r the effective 
mass is infinite! 

This same geometrical interpretation makes it possibly to 
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Figure 6 . Geometrical 
representation of the 
evolution of the natural 
forms . 

interpret the effective mass minimum of the other mode 
sometimes observed in the neighborhood of a modal inter- 
section. In practice, these properties are not observed with 
such clarity: that is because we are seldom in the presence 

of "pure" two-mode interactions and because the interactions 
with the other modes can alter this simple scheme. We will 
present in II a concrete example of such "catastrophic" 
variations of the effective mass resulting from the preceding 
interpretation . 





2. - APPLICATION TO LIQUID- PROPELLANT LAUNCHERS 

As we have indicated in the introduction, a liquid-propellant 
launcher is a typical example of an evolutive system (and 
it is within this frame of reference that we have approached 
the problem); let us recall that in the case of the Ariane 
launcher the mass of the propellants at first-stage ignition 
represents about 90% of the total mass of the vehicle! 


Moreover, we can assume with good reliability that such a 
dynamic system is symmetrical about its axis, that is to say, 
in a precise way, that it has as G symmetry group, the 
rotations of any angle about the axis of the launcher and 
all of the planes of symmetry passing through this axis. 

For a system displaying such symmetries, we know that the 
modes are classified by categories depending on the value 
of a circumferential whole number n which constitutes a 
numbering of the different irreducible representations (real) 
of 6 . In accordance with the general conclusions presented 
in I, nothing prevents the natural frequencies of modes with 
different n's to cross each other. 

In this case, we are mainly interested in the evolution of 
the vibration modes called longitudinal or axisymmetrical 
(h ~0). The corresponding irreducible reproduction of G being 
one-dimensional, these modes have a natural degeneracy of 
1. The practical interest of these modes in the discussion 
of the pogo instabilities of the launcher lies in the fact 
that they are the only ones that have a pressure fluctuation 
along the axis of the vehicle, in the immediate neighborhood 
of which the fuel lines of the propulsion systems are actually 
located. /13 

2.1. Hydroelastic vibrations: generalities 

From a mechanical standpoint a liquid-propellant launcher 
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can be represented by an elastic structure containing a liquid 
that we will assume to be incompressible and perfect. 

In order not to burden this report we will not consider here 
the effects of the stresses induced by the internal 
pressurization (reference 5) in particular, as well as the 
accoustlcal effects associated with the compressibility of 
the liquid and of the gas used for pressurization (reference 
6). In addition, we will neglect the effects due to gravity 
g and to the launcher acceleration y . 

This last approximation is fundamental to the analysis that 
will follow and therefore requires several clarifications. 
Given a hydroelastic system in an apparent gravity field 
g+y (we neglect coriolis forces); it is convenient to 
introduce two types of vibration modes defined as follows: 

The hydroelastj c modes computed by neglecting the potential 
energy terms associated with g+y. 

The sloshxng modes of the liquids computed by assuming the 
structure to be rigid and immobile. 

These two series of modes that we can consider as 
near-solutions to the entire problem must be considered as 
constituting a base for the admissible forms G of the real 
problem or, if we wish, as a system of "generalized 
coordinates" for the solution of this problem. 



The coupling between these near-solutions coming from 
satisfying the variational principle of the exact problem 
then leads to the discussion of a problem with the natural 
values of type (13). With a one-mode approach, we will say 
that the effect at the lowest order of this coupling leads 
to a renormalization of the natural frequencies of each mode; 
with a multi-mode approach, we will speak of an interaction 
between hydroelastic modes and sloshing modes. An exact pogo 
mode can then be described as a hydroelastic mode "surrounded" 
by sloshing modes. This representation of the vibration modes 
of complex structures based on the description of the resonant 
excitations in quantum physics can then be applied to the 
computation of the effects of apparent gravity Q+y on the 
pogo modes of a liquid-propellant launcher in flight or of 
modes undergoing vibration testing. 

2.2. Hydroelastic vibrations: variational formula. 

We will now show a heuristic derivation of the variational 
principle governing the hydroelastic modes in "weightless- 
ness". With regard to the rigorous derivation of the 
equations for this problem and their variational treatment, 
the reader will consult reference [8], 


To be able to use the method developed in Part I, we need 
a type (8) variational principle defined over a class of 
acceptable forms that does not depend on the fuel level of 
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the elastic structure. 


Figure 7 graphically depicts an elastic medium occupying a 
domain n in contact with a liquid with a density of pp, 
occupying a domain fip , along a surface Z. r designates the 
exposed (or free )surf ace of the liquid (we will again use 
the bold letter a to designate a vector field). 


Our objective now is to modify the Rayleigh principle (7) 
governing the vibrations of an elastic body so as to take 
into account the effect of the liquid on the wall Z of the 
elastic structure. To the denominator H 
therefore be added an additional term h 
with the displacement Up of the particles of the liquid. 


u , there must /14 

.2 


r 2 

X PF U F 


associated 


But, in the approximation of an incompressible perfect fluid, 
Up derives from a displacement potential <p (reference 8) 

Up — grad ( 15 ) 

* 

which has a modal amplitude simply related to the 
non-stationary pressure p according to (reference 8): 

p-Pp^'P (linearization of the Bernouilli Theorem) (16) 
tor g=0 

In addition, (p verifies the following equations: 

I Cup = 0 

§£-«.» over E (17) 

v 3 = 0 over r 
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expressing, respectively, the incompressibility of the liquid, 
the condition of the contact between the wall and liquid over 
Z , and the zero value of the fluctuation of the pressure 
over the exposed surface r. 

We then note that the preceding equations (17) make it 
possible to formally "express" <p "as a function" of the 
displacements of the structure along the contact surface 
Z. 

Consequently, the following equality, 

m a (u , u) = J p F m 2 „ J (grad ^2 uee (is ) 

defines a quadratic form of the displacements of the elastic 
structure. The notation has been selected to underline 
the connection between the quadratic form formally defined 
by (18) and the concept of added mass (reference 9). 




Z m Z *+ AI 


Fig. 7. -Graphical representation Fig. 8 .-Representation of 
of a hydroelastic system. the liquid domains for two 

fuel levels. 
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It is then possible to demonstrate that the stationarity 
principle 


0 = S[w 2 J = 5 



ue e 


(19) 


where [to 2 ] is a function of only the structure displacements, 
totally defines the natural forms a ^ and the natural pulses 
S associated with the hydroelastic vibration modes 

2.3. Evolution of the hydroelastic modes 

We now consider a fuel level configuration of an elastic 
structure temporarily defined by the vertical dimension of 
the horizontal exposed surface of the liquid. In the new 
configuration the liquid occupies a domain 0^, has an exposed 
surface r* and wets the structure over a surface Z*. We 
designate by means of An F ( Figure 8) the variation of the fluid 
domain bounded by the exposed surfaces r and r* and the contact 
surface AZ. 

Let M^(tt ,tt) be the quadratic form of the "added masses" for 
this new configuration. It is convenient to formally 
introduce the increase of the quadratic form between the 
configurations defined by r and r* according to: 

AA i A (u , u) - A/*(« ,u) - M a (u ,«);«€ e. ( 20 ) 
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The variational problem governing the hydroelastic vibrations 
u*, (d* of this new problem can then be written: 

2 fa V") e U (tl ^ 

0 = 5 [to* 2 ] = 6 - — > u ee, (21) 

- M a (u ,u)+ 2 Mf ' 1 C " ' M) + 2 Ja p,/2 

By looking for its u* solutions in the form of a linear 

combination Za*tt- and by noting the orthogonality relations 
j J 

(see equation (11)) verified by Uy , the "besfoty and the 
associated natural frequencies verify the following equation 
analogous to (13): 

S {/ a, oj 2 [S„ + A M a {u ( , uj)) a, = 0 ( 22 ) 

By using discretization and reduction techniques associated 
with the method of finite elements (see [8] and [6]) we can 
compute A Mi { u- , u ; ) directly and properly reach the numerical 
solution of the preceding natural value problem over a 
truncated base of natural forms u.q , Uj This method 

is being developed at ONERA for the computation of the 
hydroelastic modes of liquid-propellant launchers. It should 
make it possible to compute "by finite elements" (very 
accurate but rather burdensome method) only a limited number 
of launcher configurations and to have access to the dynamic 
characteristics for intermediate configurations through this 
simple variational method and we could also show that it is 
related to the methods called "of dynamic sub-structuring". A 6 
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Within the scope of this report, we will again simplify 
problem (19) and establish an analytical expression of-&M^, 
valid for r* configurations very close to r. 


We start from an exact expression of AM^(u,u) obtained by 
using definition (18) of and the properties of <|> and tp* 

expressed by equations (17); 


y/ ,i (« ' w ) " P p ( grad y>*) 2 - J p^Cgrad y)‘ 


~f 

t/ I 1 ♦ 


Pit* 


dp* 

bZ 


L 


(23) 


Pp'pu.n 


tp* is of course the solution of equations (17) for the same 
value of u C, but for the fl p domain bounded by T* and E*). 


We now propose not to keep the Ah terms of equation (23) and 
to neglect the higher-order terms. 


Discussion of the 
r, taking into 
r (17) , we can , 


<P with 

3^1 
32 Ir* 


a/i 

A/ 3 Z 


and 


term 


/ 
t/ M* 






bZ 


in the neighborhood of 

account the zero value condition of <p over 

within this order of approximation, replace 

neglect the difference between and 

oZ 


Discussion of the term / $ u.n : AS is "of the same order 

J AS 

as" Ah ; the same for <p ; consequently this term brings only 
negligible contributions. 

Finally, we will write for very small changes in level Ah: 


A M a (u , u) 
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(24) 



To discuss tho practical interest of this approximated 
expression, let us now consider the case of a 

liquid-propellant launcher in flight. We usually introduce 
the draining rate parameter t to characterize the fueling 
configuration of the operating stage. This x parameter is 
defined as follows: given the total mass of fuel at 

ignition of the operating stage considered, and M the mass 
of the fuels of that same stage for a certain flight 

configuration, then, by definition: 

r _ . M 
100 " A/q ’ 

(25) 

If Mj and Mg respectively designate the initial masses of 
fuels 1 and 2 (M^ = Mj + Mg ) , it will be useful, later, to 
introduce the respective mass ratios of fuels 1 and 2: 

AJ, _ A/ 2 

+M 2 ' Xj ~A/, +M 2 ' (26) 

We can now transform the analytical expression (24) of aM^ 
to approach the study of the evolution of vibrations of liquid 
launchers. By using the preceding definitions, the matrix 
elements of AM^ that enter into the solution of (22) can then 
be written, by introducing the notation d<p^}dZj to designate 
the vertical displacement of the exposed surface of the tank 
for fuel 1 for the ith mode, and <> for the average surface 
value of a quantity: /17 


OKIGI^x it.-. 
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AM A (u (l Ui)~~A T ^ r 


(27) 



str>f , 

Of* l w^t » * 



2.4 Tangent to the frequency evolution curves. 


For an infinitesimal small At «= e of the draining rate, the 
analytical expression (27) for aM^ is correct for order e , 
and we can use it within the one-mode approximation, to 
compute the frequency variation for order e , and hence, the 
derivative 6oj / Ot . By calling F the frequency w/2ji , and by 
designating with \ig the generalized mass of the mode "being 
followed" (for an arbitrary normalization), we obtains 


a f 

9 r 


JL tk 

200 M, 


(Xj <n?> + x 2 <t?5» 


(28) 


where <n?i> (resp <n 2 >.'designates the average surface quadratic 
1 2 

value of the vertical displacement of the exposed surface 
of the tank (resp 2). 


This last expression (28) can be subjected to a direct 
physical interpretation which simultaneously thrusts light 
on the approximation (24 ) (27). Indeed, we can verify that 
the terms inside parentheses give, for a variation in the 
draining rate, the amount of mass removed at the exposed 
surface of each reservoir where the modal amplitude of the 
displacements (vertical) is precisely riij and P 2 (resp). From 
this standpoint (28) is nothing but "the formula of the 
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displaced frequencies) well known to experimenters, (wo note 
the expected result oF/Ov> 0 expressing the fact that during 
the emptying of the tanks, the natural vibration frequencies 
can only increase). 

2.5. Hydroelastic modal interactions. 

Expression (27) clearly shows why two hydroelastic modes of 
the same type, let us say longitudinal for the purpose of 
illustration, necessarily "interact" during the evolution. 
Indeed, the average values of the products of vertical 
displacements of the exposed surface of each tank do not 
generally add to zero when we consider any two modes. 

2.6. Application to the flight of the Ariane launcher. 

We have numerically applied the solution to equations (22) 
simplified by the analytical approximation (27) and we have 
applied it to the case of first-stage combustion (L 140) of 
the Ariane launcher. 

Figure 9 represents the results of computations by finite 
elements (reference 10) of the asymmetrical hydroelastic modes 
of the whole launcher for four flight configurations of the 
first stage (points o). On this same figure, arrows have 
been placed to reflect the tangents to the frequency evolution 
curves computed by means of analytical expression (28) by 
using the displacements of the exposed surface and the 
generalized masses resulting from computations by finite 



elements (reference 10). 


C 


Finally, the numerical results obtained thanks to the method 
of mode interaction that we have just presented have been 
indicated by means of crosses. The continuous curves /18 
represent the most reasonable frequency evolution curves that 
we could plot considering the additional Information obtained 
by means of the modal interaction model. 



Figure 9.- Evolution diagram of the axisymmetrical vibration 
frequencies of the Ariane launcher (first-stage combustion), 
o : computation by finite elements 

+ : variational resolution of the approximate evolution 

problem. 

Considerations to be noted: 

With respect to the bold move consisting of describing 
AM by means of its linear portion (27) for such an amplitude 
in the draining rate variation, the agreement between the 
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estimates obtained from different reference configurations 
(30% and 50%) can be considered as very satisfactory, at least 
for the first hydroelastie modes. 

- We observe a "strong" modal interaction between modes 
#3 and #4 between the 50% and 70% configurations. 

Diagram 9 includes modes of a practically constant 
frequency. This means, from (28), that in practice such modes 
are not accompanied by any displacement of the exposed surface 
inside the tanks in use. The examination of the natural forms 
then makes it possible to verify that these modes are 
"localized" to the structures of the upper stages. 

These "constant" modes interact correlatively very 
little with the others (very small (u^ ,u.y ) ) : we are then 
in the presence of weak mode interactions giving rise to an 
extremely small frequency interval. Therefore, we can 
eliminate without risk such modes from the study of the 
evolution. This remark must make it possible to consider 
the computation of the pogo modes of a liquid-propellant 
launcher on the basis of pogo modes computed for a given 
configuration and selected as a function of the sufficiently 


rapid 

evolution 

criteria 

based 

on 

equation 

(28). This 

observation will 

be 

used 

profitably 

in the computations of 

finite 

evolution 

of 

the 

pogo 

modes 

by means 

of the method 

described above 

and 

which 

must 

lead 

to a 

"model with tt 


significant pogo modes" of a liquid-propellant launcher. /19 
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2.7. Evolution of the effective masses. 


(late 10 depicts the evolution of the effective mass of the 
modes (defined for a normalization along the vertical 
displacement of the first-stage engine assembly) for an 
increasing number of "test forms" injected into the 
approximate variational principle. We note, and that is 

predictable, that while taking into consideration a relatively 
restricted number of modes makes possible a reasonable 
estimate of the natural frequencies, the precise estimate 
of the evolution of the effective masses requires the 
simultaneous taking into consideration of a laige number of 
natural forms. 

These diagrams make possible the observation of effective 
mass singularities and of minimum values of for certain 
modes, in the neighborhood of a strong mode interaction. 

N.B. The evolution of the first pogo mode (at non-zero 
frequency) is sensitive to the overall translation mode at 
a zero frequency (mode #0 ) . The reader will be able to check 
that the interaction between this translation mod' i and the 
hydroelastic modes resulting from satisfying the variational 
principle (19) is directly associated with the principle of 
the conservation of momentum which, in the case 
of small motions of a free system, leads to the immobility 
of the center of gravity of this system. 
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Figure 10. - Variational extrapolation of the vibratory 
characteristics of the Ariane launcher in the neighborhood 
of configurations computed by finite elements. Evolution 
of the frequencies F and of the effective masses as a 
function of the uumber of test natural forms. 
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CONCLUSION 


During this study, we have used as theoretical cases of 
frequency crossings the case of the merging of natural 
frequencies oi modes with different types of symmetry, and 
the case of the merging of natural frequencies of certain 
modes of a system for a particular configuration corresponding 
to a sudden increase of the symmetries. 

Outside of these cases, frequency crossings can constitute 
only accidental events by comparison with the modal inter- 
actions revealed by the frequency evolution curves of the 
modes of a given type of symmetry, in the case of highly 
evolutive systems such as liquid-propellant launchers, 
missiles , etc . . . 

This study seems to us to be equally as useful for the 
analysis of the modal interactions caused by a perturbation 
brought on by the measurement itself in the case, for example, 
of the experimental technique of shifted frequencies. 

The variational method for treating evolutive conservative 
dynamic systems can be also applied to cases where the study 
of secondary effects is involved such as a change in the 
rigidity of a launcher structure resulting from tank 
pressurization (the t parameter is then the internal 
pressure ) . 
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This method could then be extended to the case of 

non-conservative systems (by dropping variational reference 
(8)), a problem of this type arising particularly in the case 
of aeroelastic instabilities of aircraft (reference 11). 

Finally, it is indispensable to mention that real structures 
do not have a perfect symmetry and that, consequently, the 
crossings of theoretical frequencies can be observed 

experimentally only as events of "weak" modal interaction 
which result, therefore, as a forbidden frequency band that 
is narrower with a smaller departure from perfect symmetries. 
These remarks must be taken into consideration to properly 
approach the problems of mode identification. / 21 
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